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We study a weakly coupled, frustrated two- leg spin- 
1/2 Heisenberg ladder. For vanishing coupling between the 
chains, elementary excitations are deconfined, gapless spin- 
1/2 objects called spinons. We investigate the fate of spinons 
for the case of a weak interchain interaction. We show that de- 
spite a drastic change in ground state, which becomes sponta- 
neously dimerised, spinons survive as elementary excitations 
but acquire a spectral gap. We furthermore determine the 
exact dynamical structure factor for several values of momen- 
tum transfer. 

PACS No: 75.10.Jm, 75.40.Gb 



I. INTRODUCTION 

The role of frustration in quasi one-dimensional mag- 
netic materials has attracted much experimental and the- 
oretical attention in recent years. On the theoretical side, 
the simplest example of a frustrated quantum magnet 
is the spin- 1/2 Heisenberg antiferromagnetic chain with 
nearest neighbour exchange SJ and next-nearest neigh- 
bour exchange J. This model is equivalent to a two-leg 
ladder (see Fig.l), where the coupling along (between) 
the legs of the ladder is equal to J (SJ). 
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FIG. 1. Heisenberg zigzag ladder. 

The zigzag ladder model is believed to describe the 
quantum magnet SrCu02 QJ^] above the magnetic or- 
dering transition, which takes place at about T ss IK. 
The exchange constants are estimated to be J w 1800.fr, 
\5J/J\ ~ 0.1 — 0.2 Q. A second material with zigzag 
structure, that has recently attracted much interest, is 
CS2CUCI4 Hi . However, in CS2CUCI4 all neighbouring 
chains are coupled by a zigzag interaction and no pro- 
nounced ladder structure exists. 



In Refs. jl|-f§ it was argued that a weak antiferromag- 
netic zigzag coupling between the chains drives the model 
to a massive phase, characterised by spontaneous dimeri- 
sation (see also [Q). Let us briefly review some important 
parts of the derivations of |^,^| . The lattice Hamiltonian 
of the zigzag ladder is 



j=l,2 n 

+ SJ (Sl,n + Si, n _|_i) • S2 in , 



(1) 



where we assume that 5 -C 1. The low-energy effective 
action for ([j]) is now obtained as follows. For 6 — ► one 
is dealing with two decoupled Heisenberg chains, which 
can be bosonised in terms of two Wess-Zumino-Novikov- 
Witten (WZNW) models by using the standard relation 
between the spin density on chain j and the fields of the 
WZNW model (see e.g. f§) 



[jf{ X )+jf(x)]+{-ir^(x) 



(2) 



Here ao is the lattice spacing, and the fields J" and J" 
are the right and left currents of the WZWN model cor- 
responding to chain j. They parametrise the smooth 
component of the magnetisation. Finally, rij is the stag- 
gered component of the magnetisation on chain j. Using 
(||) , the zigzag interchain interaction can be expressed in 
terms of the WZNW fields. In this way one straightfor- 
wardly obtains the current-current interaction J5|,|| 



W c = Ai(J f + J?)(J 2 ° + Jf) 
— Ao(^i J\ + J% J2 ) ' 



(3) 



where Ai cx SJ. A standard renormalisation group (RG) 
analysis then shows that the antiferromagnetic interchain 
interaction Ai leads to a spontaneously dimerised ground 
state In |j] it was shown that, in addition to the 

current-current interaction (^), a "twist" term arises 



(4) 



In the presence of exchange anisotropics the twist term 
induces incommensurabilities in the spin correlations |j) . 
We expect this to be hold true even in the SU(2) sym- 
metric case (no exchange anisotropics) we are interested 
in here. In the latter case it can be shown that the twist 
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term and current-current interaction are equally impor- 
tant in the RG sense: they diverge (i.e. reach strong 
coupling) simultaneously, with a fixed ratio pC[ . As far 
as the SU(2) symmetric zigzag ladder is concerned, it is 
therefore not possible to separate the effects of the twist 
and current-current interactions in a simple way. 

However, from a purely theoretical point of view it 
clearly is desirable to develop a thorough understanding 
of the physics due to isolated current-current and twist 
interactions. Their effects can be disentangled by intro- 
ducing an exchange anisotropy || , which makes the twist 
more and the current-current interaction less relevant in 
the RG sense. Using this trick, a pure twist interaction 
was studied in ||. 

The role of an isolated current-current interaction 
has been previously investigated in connection with the 
zigzag ladder in [0J|^]. In particular, the spectrum of 
elementary excitations and the dynamical structure fac- 
tor were calculated in |(| using large-N techniques. It is 
known, that extrapolation of large-N results to small val- 
ues of N can lead to incorrect results Jll]] . Having this in 
mind, we carry out an exact calculation in order to de- 
termine the spectrum and structure factor. We find that 
the large-N results are indeed qualitatively incorrect. 

As explained above, in the zigzag ladder both twist 
and current-current interactions are present; therefore, 
strictly speaking, our results cannot be directly applied 
to this model. Nevertheless, we believe that many of 
our findings presented below remain qualitatively correct 
when applied to model (0) . We discus this point in more 
detail in section VIl . 



In order to connect our results to a microscopic model, 
we consider a frustrated spin ladder modified in such a 
way that only current-current interactions emerge in the 
low-energy effective action. 

The outline of this paper is as follows: in section || we 
introduce a frustrated spin ladder model giving rise to the 
desired low-energy effective field theory. In section |n| we 
show that the resulting field theory is essentially equiv- 
alent to an 0(4) Gross-Neveu model Sections IV 
and are concerned with the description of the ground 
state(s) and elementary excitations. In section VI we de- 
termine the exact dynamical structure factor for several 
values of momentum transfer and show that there are 
no coherent contributions to the structure factor. We 
conclude with a summary and discussion of our results. 



II. A FRUSTRATED LADDER WITHOUT TWIST 

The model we consider is a generalisation |l3| of the 
standard two-leg spin ladder which, apart from the on- 
rung coupling J± , also includes an interaction J x across 
both diagonals of the plaquettes. The Hamiltonian reads 



H — J Sj >n ■ Sj,n+l + Jj_ Si >n ■ S 



2.0 



] (5) 

n 

We assume that 

J, J±, J x > 0, J> J x , Jx • (6) 
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FIG. 2. The twistless ladder model. 

The low-energy effective action can be derived by non- 
abelian bosonisation in the usual way. The Hamiltonian 
density is found to be of the form 



H(x) =Hi(x) + H 2 (x) +Hint[x) , 



where TL1.2 are critical SU\{2) WZNW models with a 
marginally irrelevant current-current perturbation (Ao > 
0): 



n, 



2lTV* 



3 ^--j -> 3 
-A Jj • Jj , j = 1,2. 



J, -J, 



(8) 



The interaction part is given by 

Hint = Ai (Ji + Ji) • (J 2 + J 2 ) + A 2 n x • n 2 , (9) 
with the coupling constants 



Ai = (Jj_ + 2J x )a , A 2 = (Ji - 2J x )a . 



(10) 



No marginal perturbation with the twist-term structure 
arises because the staggered magnetisation operators add 
rather than subtract due to the geometry of the problem. 
The absence of such term can also be deduced from the 
existence of discrete (reflection) symmetries of the lat- 
tice Hamiltonian If J± = 2J X only the marginal 
(current-current) interaction survives. This is the case 
we study in the remainder of this paper. 

We note that for generic values of J± and J x the in- 
teraction of staggered magnetisations dominates and the 
resulting physics is essentially the same as for the stan- 
dard ladder (J x = 0) @ (see also chapter 21 of @). 



III. DUALITY TRANSFORMATION 

The low-energy effective action (||)-(^|) can be recast 
as a theory of four massive, interacting, real (Majorana) 
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fermions, or equivalently, four weakly coupled Ising mod- 
els§ 

. 3 

H = - ^ V a (4>ad x 4> a - 4> a d x 4>a) 
a=0 

X± — A 3 _ _ Ai + A - . 3 „ 
H 2 E ^^j^j o V'oV'o^V'tV'i- ( n ) 

j>i— 1 i—l 

Here v± — v 2 = V3 = v s ^ vq are the velocities of the 
four Majorana fermions. The lattice spin operators are 
expressed in terms of the Majorana fields and order and 
disorder operators of the four Ising models as 

S z + (x) cx -i (V-i^ + &$a) ~ -4(-l) x/o >iM2^3^o , 

S z _(x) <x i (tfsVo + Mo) + A{-l) x/ao a l( j 2 ^ , (12) 

where S±(x) = Sf(x) ± S 2 (%) and A is a nonuniver- 
sal constant. Analogous expressions are available for the 
other components of the spin operators [^) . A standard 
one-loop RG analysis shows that the coupling Ao flows 
to zero, so we will ignore it in what follows. In order to 
further simplify the problem, we also neglect the small 
difference between the velocities v s and vo, and finally 
perform a duality transformation on the O-Majorana 

ipo -> ip4, , ^0 -> ~i>i , ct -> m , Mo -> cr 4 . (13) 

This yields the Hamiltonian of the 0(4) Gross-Neveu 
model (H 

4 A 4 

W = — ^ ^id x 1pi - 4>id x 1pi + — ^Pi^i^j^j ■ 

i=l 3>i=l 

(14) 

Under @) the spin-densities transform to 

(x) cx -i (^1^2 + M2) - ^(-l) x/oo MiM2<T3M4 , 
S* (1) (x z (^4 - ^3^4) + A{-l) x/aa a 1( j 2 ^a A . (15) 



A. Twistless Ladder 

The low-energy effective model ( |l6| ) exhibits a local 
Z 2 symmetry related to independent translat ions by one 
lattice spacing on each chain [tp± —* tp± + yJir/2). This 
symmetry is spontaneously broken in the ground state 
and leads to a nonvanishing dimerization. Notice that 
the Z 2 symmetry appears to be a feature of the low- 
energy sector only and follows from the fact that spin 
currents J1.2 are translationally invariant objects. The 
transformation Si(n) — > \ [Si jn +i + Si )n _i], or a similar 
one with Si^. — > S 2 ,k, changes the lattice Hamiltonian 
but leaves the low-energy effective field theory invariant 
and maps the two ground states onto one another. 

In order to characterise the dimerisation patterns of 
the two ground states, we determine the expectation val- 
ues 

(Si, n ■ S 2 , n ) oc (Ji{x) ■ J 2 (x)) + (ni(x) ■ n 2 (x)) , 
(Si,n ■ S 2 ,n+i) oc (Ji(x) ■ %i{x)) - (ni{x) ■ n 2 (x)) , 
(5*2, n ■ Si.n+i) oc (Ji(x) ■ J 2 (x)) - (n\{x) ■ n 2 (x)) . (17) 

After performing the duality transformation to the 0(4) 
Gross-Neveu model and bosonising, we obtain 

(fti(x) ■ 7x2(2;)) cx (cos v / 27rcp+ cos \p2/Kkp-) = ±const m , 




= const m 2 , (18) 

where m oc exp(— const J/ J±) is the (exponentially small) 
soliton mass in the Sine-Gordon model. Due to the small- 
ness of m, the (n\(x) ■ n 2 {x)) expectation value domi- 
nates in (|l7j), so that within the exponential accuracy 
the dimerisation is proportional to the quantum soliton 
mass. 



IV. GROUND STATE 



In order to proceed, it is convenient to use the rep- 
resentation of ( |l4] ) in terms of two sine-Gordon models 
[pTsfl . Ignoring terms that only renormalise the velocity 
we find that ( |i"4| ) is equivalent to 



i=± 



+2Ai 



(27ra ) 2 



(16) 



where Qi are the dual fields. The two sine-Gordon models 
( |l6| ) occur on the SU(2) invariant strong-coupling sepa- 
ratrix of the Kosterlitz-Thouless phase diagram and are 
thus in the massive regime. 




FIG. 3. Qualitative picture of the two degenerate dimerised 
ground states: spins connected by the solid (dotted) lines have 
a tendency to form singlets (triplets). 
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The Z 2 symmetry of the low-energy effective Hamilto- 
nian ([If]), that manifests itself in the degeneracy of the 
two ground states corresponding to different signs in (18), 
is spontaneously broken, implying the existence of mas- 
sive Z 2 kinks. It turns out (see below) that these kinks 
are elementary excitations of the model. 



B. Zig-Zag Ladder 

Let us discuss the implications of the emergence of 
spontaneous dimerisation for the case of the zig-zag lad- 
der if we ignore the twist term. For the zig-zag ladder 
the appropriate definition for the dimerisation is 

d=ls x {x)- (s 2 (x + ao/2)-S 2 (x-a /2))y (19) 

In the continuum limit we find 

d oc (cos cos ^/2~nipJ) — zbconst m . (20) 

The resulting dimerisation patterns are shown in Fig.^J 

WWW 
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FIG. 4. Qualitative picture of the spin configuration in the 
dimerised ground states: spins along the thick diagonal bonds 
have a tendency to form singlets. 

We believe that taking into account the twist term will 
not qualitatively change this picture. 



V. EXCITATIONS 

From the exact solution of the Sine-Gordon models 
( |l6| ) we infer that there are only four elementary exci- 
tations corresponding to solitons and antisolitons in the 
± sectors. We denote these by s± and s±. The ele- 
mentary excitation have a simple interpretation in terms 
of dimerisation kinks, i.e. domain walls separating re- 
gions of dimerisation with opposite sign. It can be shown 
along the lines of Q that these particles carry spin ±1/2. 
In terms of the low-energy effective theory of two Sine- 
Gordon models (16), the total spin density is given by 



Sl{x) + S%(x) = -L [dztp+ix) + d x <p-(x)] . (21) 



/2tt 

Kinks interpolate between asymptotic values of the fields 
ifi differing by ± y/n/2 as is most easily deduced from the 
fact that the classical vacua of (llq) are located at 



i I class 



± 



(22) 



where are arbitrary integers. Integration of (21) then 
yields that a single kink carries spin 



± 



1 



2tt V 2 2 



(23) 



The results presented below for the dynamical struc- 
ture factor are consistent with the interpretation of these 
kinks as gapped spinons. Altogether there are two spin- 
1/2 multiplets, corresponding to one multiplet for each 
leg of the ladder. The emerging physical picture is quite 
simple and pretty: the two-spinon states observed in the 
structure factor simply correspond to the kinks related to 
the spontaneous breakdown of the discrete Z 2 symmetry. 

Simple visualizations of this picture are shown in Fig.|5| 
for the twistless ladder and in Fig.fj] for the zig-zag ladder. 




FIG. 5. A two-spinon state in the twistless ladder. Spinons 
correspond to kinks between domains with different sign of the 
dimerisation. Solid lines depict bonds along which there is a 
tendency to form singlets. 

For the twistless ladder the kinks correspond to verti- 
cal domain walls between regions with different signs of 
dimerisation. There is a spin- 1/2 associated with each 
domain wall, although this is not immediately obvious 
from Fig.^J. In order to get a feeling why a spin- 1/2 might 
be associated with each kink, let us think of the trans- 
lationally invariant, "double-zigzag" ground state shown 
in FigJ^ as a symmetric superposition of two dimerised 
states. Each such state represents a sequence of plaque- 
ttes with ideal singlet bonds across the plaquette diag- 
onals (with each spin involved in one bond only), has a 
period 2ao and is shifted with respect to the other state 
by one lattice spacing. If the "double-zigzag" phase occu- 
pies a finite domain of the ladder, for the two 2ao-periodic 
dimerised states to resonate, the number of rungs within 
such a domain should be odd. Then the two-kink config- 
uration in Fig.^ can equivalently be viewed as the super- 
position of states shown in Fig.|[ 



4 




S(u>, q, q±) oc Imi 



dx / dte~ iut+iqx 
) Jo 

x([{S* 1 (t,x)±SI(t,x)},{S* 1 (0,0)±SI(0,0)}}) • (24) 



where the positive (negative) sign corresponds to q± = 
(<?_L = tt). 



A. Summary of large- N results 




FIG. 6. "Resonating" ideal dimer configurations. 

The intuitive picture one obtains from Fig|j| is then 
that on average there is indeed a spin- 1/2 associated with 
each kink. 

For the zig-zag case a much nicer picture emerges. The 
Z2 symmetry corresponds to a reflection symmetry on the 
lattice and kinks look like left over spin-l/2's as shown 
in Fig.g 

W/AV 

FIG. 7. Physical picture of a two-spinon state. Spinons 
correspond to kinks connecting domains with different sign of 
the dimerisation. 

The intuitive picture of Fig. Q fits well to the identi- 
fication of a spinon in a spin- 1/2 chain as a bare spin 
insertion into the ground state |17f| . 



The dynamical structure factor has been previousl 
calculated in the framework of a large- N approach | 
The limit N — > 00 of (jl4|) is equivalent to a theory of free 
massive Majorana fermions 



N 



(25) 



The presence of the mass term reflects the spontaneous 
breakdown of parity, which in turn implies the existence 
of two degenerate ground states. The sign of the mass 
terms, as well as the expectation values (a a ) and (/i a ), 
depend on the choice of ground state. In the case where 
{a a ) 7^ (the mass of the triplet is positive), the structure 
factor for q± = and q around 0, it was shown to be || 



S(u, q ss it, 0) oc - — -5 (u — \/v 2 s (q~ it) 2 + m 2 ) , 



S(u,q*i 0,0) oc 



2 2 
771 q 



s 3 Vs 2 — 4m 2 



(26) 



where s 2 



or 



2 2 
%q- 



The explicit expressions for the 



structure factor around (q,q±) — (0,7r), (it, it) are com- 
plicated, but reveal the presence of incoherent two and 
three particle continua, respectively. We will now show 
that the results obtained in the large- N limit are qualita- 
tively incorrect. The reason for this failure of the large- 
N approach is that it entirely neglects the existence of 
topological kinks interpolating between the two degener- 
ate ordered ground states. Extrapolation of the large- iV 
results to lower values of N shoud be done with caution 
because the spectrum of the O(N) Gross-Neveu model is 
very sensitive to the value of N DTl . 



VI. DYNAMICAL STRUCTURE FACTOR 

The long-distance asymptotics of the spin-spin cor- 
relation functions are dominated by the soft modes at 
q = 0,7r, q± = 0, 7r, where q and q± denote the wave- 
numbers along and perpendicular to the two chains, re- 
spectively. In what follows we will determine the dy- 
namical structure factor for wave numbers in the vicinity 
of the above four points in g-space. Due to the spin- 
rotational symmetry the dynamical structure factor is 
given by 



B. Exact results 

We now determine the dynamical structure factor us- 
ing exact results on formfactors in the Sine-Gordon model 

jnj-§II- We start with the case 9-l = °. 1 ~ °- The 

smooth component of the sum of the two spin densities 
is expressed in terms of the Sine-Gordon models as fol- 
lows 



Si(x)+SZ(x) 



(x d x ip+ + d x ip. 



(27) 



oth 
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This is nothing but the sum of the temporal components 
of the current operators in the two Sine-Gordon models 
(j+ +3-)- We are interested in the structure factor, i.e., 



S(u>, q ps 0, 0) oc Im i / dx 



fl£ ^i(ui+ie)t — ivqx 



x(\j° a (x,t) , j°(0,0)]> 



(28) 



where v s is the velocity of the excitations. We express 
( p8| ) in the spectral representation using our knowledge 
of a complete set of states in terms of (anti) soliton scat- 
tering states. Energy and momentum are parametrised 
in terms of the rapidity variable as 



p = in sinh 9 e = m cosh 9 , 



(29) 



where m is the mass of the four elementary excitations. 
The resolution of the identity is given by 



n—0 a>i 



dB-L ■ ■ ■ d9 r , 
(27r)"n! 



h)a n --- ai ai ' a " (#1 ' 



(30) 



where n is the number of particles and on € {s±,s±} 
specifies their respective "flavour" (soliton or antisoli- 
ton in + or — sector). Inserting ( |30| ) in ( p8|) and using 
Poincare invariance yields 

S{u,qK 0,0) tx 

n=0 at \ ) ■ 

5 ( v s q — tci sinh 9 A 5 (v s q + m sinh 9j 



— m V ■ c °sh 9j + ie 



lo + m J2j cosh 9j + ie 



(31) 

where Fjo(9i ■ ■ ■ 9 n ) ai ... an is the Sine-Gordon current 
form factor 



F j o(9 1 ---9 n ) ai ... an = (O|j°(O,O)|0„ 



7 l/a n ■ --CKi • 



(32) 



We note that an n-particles state only contributes to (pi]) 
above the n-particle threshold, i.e. s 2 — uj 2 — v^q 2 ^^ 
n 2 m 2 . Thus, at low energies s 2 < 16m 2 only two-particle 
states contribute. The corresponding formfactor is 



F j0 (9 1 ,6 2 ) S s = -2m sinh ( ) f(9 1 



.sinh 6»/2 




, /7TK\ 



K sh(7TK) 

After performing the ^-integrations we obtain 



(33) 



S(lo, q re 0,0) cx 



2 2 2 



l/(20( S ))| S 



s 3 -\/s 2 — 4m 2 



(34) 



where 0(s) = arccosh (if-) and 4m 2 < s 2 < 16m 2 . As 
we already mentioned, the result (|34|) is exact as long as 
s 2 < 16m 2 . For larger energy transfers there are (small) 
corrections due to four, six, eight etc particle states. 
These can be calculated in the same way as the two- 
particle contribution. Approaching the threshold s = 2m 
from above, (|34|) goes to zero like y/ s — 2m. 

The result ( |34[ ) has the same structure as the one ob- 
tained in the large-A approximation. We note that the 
vanishing of the structure factor for q = (S(u>, 0, 0) = 0) 
reflects the fact that the z-component of spin is a con- 
served quantity. 

Next, we consider the structure factor at (q ss 0, 7r). 
The smooth component of the difference of spin densities 
is 



oc d T ip+ — d T ip- 



(35) 



smooth 



This is precisely the difference of the spatial components 
of the currents in the two Sine-Gordon models — ji ) . 
Using the exact two-particle formfactor we obtain the 
leading contribution to the structure factor 



S(u, q pa 0, it) oc 



m 2 m 2 |/(20( S ))| 2 
s 3 V s 2 — 4m 2 



(36) 



where f(6) is given by ( |33| ) and again 4m 2 < s 2 < 16m 2 . 
Note that the structure factor does not vanish for q — > 
as the magnetisation difference between chains is not con- 
served. This is due to the fact that our starting point 
does not have 0(4) symmetry: after the duality trans- 
formation we obtain an 0(4) symmetric Lagrangian, but 
correlation functions transform nontrivially. This result 
is of course expected, since the interchain interaction 
must break the 0(4) - SU(2) x St/ (2) down to SU(2). 

Finally, we examine the structure factor at (g ps 7r, 0) 
and (q ps it, it). The bosonised forms for the staggered 
components of the sum and difference of the spin densi- 
ties arc found to be 



Sf(x) + SI(x) 
Sf(x) - S!(x) 



stagg 



M.agg 



cos a/tt^ cos VttQ, 



sin v / tt < I ) sin ^/ttQ , 



(37) 



where $ = (ip+ + tp-)/V% and 6 = (0+ - 0_)/\/2. At 
present it is not known how to calculate formfactors for 
the operators appearing in ( ^7[ ) as they involve both the 
field and the dual field. However, it is still possible to de- 
termine the qualitative behaviour of the structure factor. 
From (|37]) it is clear that the structure factor involves the 
calculation of formfactors of operators 



cos or sin 



cos or sin 



6 



x cos or sin 



cos or sin 



(38) 



These formfactors are obviously products of formfactors 
in the two Sine-Gordon models. Let us therefore concen- 
trate on the + sector for the time being. It was shown 
in |^2j that the operators cos \f\Q+ and sin yf^0+ in 
the Sine-Gordon model with coupling constant /3 = y/8n 
have fermionic character and thus have nontrivial form- 
factors with one-soliton states. On the other hand, we 
know from that cos y/~^ip + and sin ^J^ip + are of 
bosonic character. We therefore conclude that + part of 
the operator (|3^) has fermionic character. This implies 
that it couples only to states with at least one (anti) 
soliton. An analogous statement holds true for the — 
sector, so that the leading contribution to the structure 
factor comes from two-particle states. In other words no 
coherent one-particle excitation exists. 

From the above results for the dynamical struture fac- 
tor we deduce that the low-lying excitations are described 
in terms of a gapped two-particle scattering continuum. 
As we have mentioned above, the elementary excitations 
carry spin-1/2. This leads us to identify them as massive 
spinons. 



VII. SUMMARY 

We have studied the effects of pure current-current in- 
teractions in a frustrated two-leg spin ladder. We have 
shown that spinons, which are gapless topological exci- 
tations propagating along decoupled Heisenberg chains, 
survive as elementary excitations in the frustrated ladder, 
but acquire a finite mass gap. We have given an inter- 
pretation of these massive spinons as quantum dimerisa- 
tion kinks. The kinks are deconfined and, in all physical 
states, appear only in pairs. As a result their contribution 
to the dynamical structure factor is entirely incoherent. 
Our findings bear a strong resemblance to those of |Q . 

We believe that our results not only apply to the lad- 
der (Eh, but with some modifications also to the zigzag 
ladder ([!]). As discussed above, in the zigzag case there 
is a twist term in addition to the current-current inter- 
action. We conjecture, that the effect of the twist term 
is merely to shift the minimum of the two-spinon con- 
tinua at (q = tt, 0) and (q = n, tt) to incommensurate 
wave numbers, i.e. to (q = tt + 8, 0) and (q = tt + 5, tt) , 
where \5\ <C 1. Such a picture is consistent with what 
is known from numerical studies |5|,p^-p6| and also fits 
well to what one would expect on the basis of an (uncon- 
trolled) extrapolation of the results for 6 = 0(1) [ p7|j28| 
to \S\ < 1. 

Coming back to the twistless chain (|5|), it should be 
pointed out that its ground state and excitations have 
been previously studied for the special case J x = J 
p9| , [30f ("Bose-Gayen model"). In this case, the Hamilto- 
nian (|5|) exhibits an enlarged (local) symmetry, related to 
the interchange Si(n) <-» S2(n) at arbitrary rung n, and 



decouples into two commuting parts describing either an 
array of entirely decoupled on- rung singlets or an effective 
S=l chain Q . In both cases, the ground state belongs to 
the universality class of the (undimcrised) Haldane spin 
liquids with the spin-1 massive magnons being coherent 
elementary excitations jl4| , |23]| . This is in marked con- 
trast with our fmidings for J x = |Ji « J and implies 
the existence of a crossover between the two regimes at 
some intermediate coupling. 

It should be understood that the region where the 
marginally perturbed ladder (A 2 = 0) and the Bose- 
Gayen model start overlapping, i.e. the vicinity of the 
point J_l = 2J X = 2 J, is not accessible within our contin- 
uum approach, based on the assumption that J±,J X <C 
J. Staying on the line J± = 2 J x and increasing J x would 
enforce the amplitude of the current-current perturbation 
(Ai) to increase, in which case no reliable conclusions are 
available. On the other hand, one can start approaching 
the Bose-Gayen regime by keeping Jj_ fixed and increas- 
ing J y . In this case one inevitably deviates from the 
line Jj_ = 2 J x , and that gives rise to the appearance of 
the strongly relevant perturbation A2II1 • n 2 . The latter 
introduces an extra potential, 



U 



A2[2C0S V27T (<p + — (f_) — COS V27T (if + + tpj) 

+ cosV2tF(0 + -0_)], (39) 



that couples the two Sine-Gordon models (jig), removes 
the Z2 degeneracy between the two dimerised ground 
states and thus leads to soliton confinement. The soliton- 
antisoliton pairs start forming triplet and singlet massive 
bound states and transform to coherent single-particle 
excitations. If the deviation from the line Jj_ — 2J X 
is large enough, the A2 perturbation takes over, and the 
effective low-energy field theory becomes that of four Ma- 
jorana fermions, with a mass term 



iA 2 ^2 i>o.^a - 3?Poipo 



as it is the case for the standard (nonfrustrated) lad- 
der 0. It is therefore tempting to speculate that the 
two massive Haldane phases on the both sides of the line 
J x = 2J x can be smoothly connected with those of the 
Bose-Gayen model. This, however, does not exclude the 
existence of other phases in the 3-parameter space of the 
model (||). 

As discussed in Refs. ]3l| , p8| ], a similar soliton confine- 
ment scenario is realized if one adds an explicit dimeri- 
sation to the zigzag Hamiltonian (|l]). 
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